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This  document  presents  both  the  written  and  oral  versions  of  a presenta- 
tion given  at  the  1977  IEEE  International  Conference  on  Acoustics,  Speech,  and 
Signal  Processing,  in  Hartford,  Connecticut,  10  May  1977.  The  presentation 
concerned  optimum  element  placement  for  passive  localization. 

For  mathematically  tractable  assumptions,  arrays  with  a third  of  their 
limited  elements  at  each  end  and  the  middle  of  a line  are  in  some  sense  best— 
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for  simultaneous  range  and  bearing  estimation, 
spaced  at  half  wavelengths. 
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INTRODUCTION 

The  purpose  of  this  talk  is  to  discuss  the  passive  technique 

OF  ESTIMATING  THE  RANGE^  AND  BEARING^  By  TO  AN  ACOUSTIC  SOURCE 
USING  M HYDROPHONES  IN  A COLLINEAR  ARRAY  OF  LENGTH^  L.  I WILL 
DISCUSS  THE  MAXIMUM  LIKELIHOOD  ESTIMATE  FOR  RANGE  AND  BEARING^  THE 
VARIANCE  OF  THE  RANGE  AND  BEARING  ESTIMATES^  AND  THE  OPTIMUM  HYDRO- 
PHONE PLACEMENT  FOR  ESTIMATING  BOTH  RANGE  AND  BEARING.  An  IMPORTANT 
TERM  IN  THESE  RESULTS  WILL  BE  THE  EFFECTIVE  ARRAY  LENGTH^  L SUB  Ey 
WHICH  IS  EQUAL  TO  L TIMES  SINE  B. 

- SLIDE  ONE  PLEASE  - 
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The  particular  geometry  of  interest  is  two  dimensional  with 

AN  ACOUSTIC  POINT  SOURCE  WHOSE  RANGE  AND  BEARING  ARE  TO  BE  ESTIMATED 
BY  A FIXED  NUMBER  OF  RECEIVERS.  FoR  THE  PURPOSES  OF  THIS  WORK>  WE 
PRESUME  THAT  THE  RECEIVING  HYDROPHONES  ARE  COLLINEAR.  HOWEVER, 
REGARDLESS  OF  THE  HYDROPHONE  POSITIONS,  A FIXED  NUMBER  OF  SENSORS 
HAVE  AN  INHERENT  UNCERTAINTY  IN  ESTIMATING  SOURCE  LOCATION.  ThIS 
UNCERTAINTY  REGION  IS  NOMINALLY  ELL I PT I CAL; SO  THAT  BY  PROPERLY 
DEFINING  HOW  RANGE  AND  BEARING  ARE  MEASURED.  THE  ESTIMATION  ERRORS 
CAN  BE  DECOUPLED.  FOR  A COLLINEAR  ARRAY  OF  SENSORS, WE  MEASURE  THE 
BEARING  AS  THE  ANGLE  BETWEEN  THE  LINE  ARRAY  AND  THE  MAJOR  AXIS  OF 
THE  UNCERTAINTY  REGION.  An  IMPORTANT  CONSIDERATION  IN  DETERMINING 
THE  SOURCE  LOCATION  IS  THE  EFFECTIVE  LENGTH  OF  THE  ARRAY,  L , AS 
SEEN  BY  THE  SOURCE.  FOR  THE  PURPOSES  OF  OUR  ANALYSIS,  THE  SOURCE 
IS  ASSUMED  TO  BE  SPATIALLY  STATIONARY,  THAT  IS,  NOT  MOVING  OVER  THE 
OBSERVATION  PERIOD.  ThE  TALK  BY  PROFESSOR  KNAPP  THAT  FOLLOWS  THIS 
DISCUSSION  DEALS  WITH  SOME  OF  THE  PROBLEMS  PRESENTED  BY  MOVING 

SOURCES.  Implicit  in  the  discussion  that  follows  are  Gaussian 

SIGNALS  AND  NOISES  THAT  ARE  MUTUALLY  UNCORRELATED.  RELATED  WORK 
IS  TREATED  BY  McDoNALD  AND  SCHULTHEISS,  HaHN,  AND  BaNGS  AND 
SCHULTHEISS. 


- NEXT  SLIDE  PLEASE  - 
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For  a specified  array  geometry^  the  maximum  likelihood 

ESTIMATE  OF  RANGE  AND  BEARING  IS  OBTAINED  BY  COHERENTLY  PROCESSING 
THE  OUTPUTS  OF  THE  SENSING  HYDROPHONES  AS  SHOWN  ON  THIS  DIAGRAM. 

In  particular^  each  hydrophone  output  is  prefiltered  to  accentuate 

HIGH  SIGNAL-TO-NOISE  RATIO,  THEN  DELAYED  AND  SUMMED.  ThE  SUMMED 
SIGNAL  IS  FED  TO  A FILTER,  THEN  SQUARED  AND  AVERAGED  FOR  THE 
OBSERVATION  TIME.  ThE  OUTPUT  OF  THIS  NETWORK  IS  MAXIMIZED  THROUGH 
E INDIRECT  ADJUSTMENT  OF  THE  DELAY  PARAMETERS.  ThE  DELAY  PARAM- 
,i.RS  ARE  DERIVED  ON  THE  BASIS  OF  TWO  ADJUSTABLE  PARAMETERS 
HYPOTHESIZED  BEARING,  B TILDA,  AND  HYPOTHESIZED  RaNGE,  R TILDA. 

Thus,  an  operator  need  only  adjust  his  best  estimate  of  bearing 

AND  RANGE.  FrOM  THESE  TWO  INPUTS,  PROPER  DELAYS  ARE  INSERTED  IN 
EACH  HYDROPHONE  RECEIVING  LINE.  ThE  PROCESS  OF  DELAYING  AND 
SUMMING  IS  A GENERALIZED  BEAMFORMER  WHERE  THE  DELAYS  USED  CAUSE 
THE  BEAMFORMER  TO  PRESUME  THE  SOURCE  WAVEFRONT  IS  CURVED  AND  NOT 
PLANAR.  The  individual  SENSOR-TO-SENSOR  delays  inserted  are  DI- 
RECTLY RELATED  TO  THE  HYPOTHESIZED  SOURCE  AND  SENSOR  LOCATIONS. 

-NEXT  SLIDE  PLEASE- 
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There  are  several  physical  interpretations  for  the  maximum 

LIKELIHOOD  ESTIMATES  FOR  RANGE  AND  BEARING.  OnE  INTERPRETATION 
IS  SHOWN  HERE.  By  HYPOTHESIZING  A SOURCE  POSITION  CHARACTERIZED 
BY  A FINITE  RaNGE,  R TILDA^  AND  BEARING,  B TILDA,  A SPHERICAL  OR 
CURVED  WAVEFRONT  IS  PRESUMED.  ThIS  CURVATURE  AND  THE  ACTUAL 
HYDROPHONE  POSITIONS  GIVE  RISE  TO  DELAYS  THAT  ARE  INSERTED  IN  THE 
GENERALIZED  BEAMFORMER  TO  MATCH  THE  PROCESSOR  TO  THE  WAVEFRONTS 
CURVATURE. 

- NEXT  SLIDE  PLEASE  - 
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Cramer  Rao  Lower  Bound  (Min.  Var.) 


27tC^ 

T W((jj)<Ai^6<jj 


Where 

T is  observation  time 

is  a function  of  the  sensors  and  their 
positions  (Different  for  R,  B) 

V is  a function  of  the  output  SNR 

C is  the  speed  of  sound  in  the  medium 
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SLIDE 


Bangs  and  Schultheiss  have  shown  that  the  Cramer  Rao  lower 
BOUND  for  estimating  A PARAMETER  « IS  CHARACTERIZED  BY  THE  SIMPLE 
FUNCTION  SHOWN  HERE  (WHEN  S I GNAL-TO~NO I SE  RATIO  AT  EACH  SENSOR  IS 
THE  same).  The  parameter  $ can  denote  range  OR  BEARING.  In  par- 


ticular^ THE  VARIANCE  OF  THE  ESTIMATED  PARAMETER  IS  CHARACTERIZED 


BY  TWO  PI  TIMES  C SQUARED  WHERE  C IS  THE  SPEED  OF  SOUND  IN  WATER; 

OVER  THE  PRODUCT  OF  THREE  TERMS.  ThE  THREE  DENOMINATOR  TERMS 
INVOLVE  Uj  AND  V TO  BE  DISCUSSED  ON  THIS  AND  SUBSEQUENT  SLIDES:  ) 

I IS  THE  OBSERVATION  TIME.  | 

U SUB  £ IS  A FUNCTION  OF  THE  SENSORS  AND  THEIR  POSITIONS. 

For  variance  of  the  range  estimates  we  use  U sub  R and  for  j 

VARIANCE  OF  THE  BARING  ESTIMATES  WE  USE  U SUB  B. 

Y IS  A FUNCTION  OF  THE  S I GNAL-TO-NO I SE  RATIO  AND  THE  NUMBER 
OF  SENSORS.  Note  that  it  is  the  integral  of  V of  omega  over 

THE  FREQUENCY  BAND  ZERO  TO  CAP  OMEGA  THAT  IS  IMPORTANT  IN 
THE  VARIANCE  EXPRESSION.  AlSO^  HIGH  FREQUENCIES  ARE  MORE 
IMPORTANT  THAN  LOW  FREQUENCIES,  DUE  TO  THE  OMEGA  SQUARED 
TERM  IN  THE  INTEGRAL.  ThE  CrAMER  RaO  BOUND  PRESENTED  HERE 
IS  FOR  LOCAL  VARIATION  AND  IMPLICITLY  ASSUMES  HIGH  OUTPUT 
SIGNAL-TO-NOISE  RATIO. 


Increases  in  the  product  of  T,  U,  and  V result  in  decreases  in 

THE  VARIANCE  OF  THE  RANGE  OR  BEARING  ESTIMATES.  FOR  THE  ARRAY 
CONFIGURATIONS,  T,  U,  AND  V ARE  INDEPENDENT;  SO  INCREASING  T,  U, 
OR  V WILL  RESULT  IN  DECREASED  VARIANCE.  FOR  EXAMPLE,  DOUBLING  T 


■I 
■ > 

1 

4 


WILL  REDUCE  THE  STANDARD  DEVIATION  OF  EITHER  THE  BEARING  ESTIMATE 
OR  RANGE  ESTIMATE  BY  ABOUT  ONE  POINT  FOUR.  MINIMIZE  THE  VARIANCE 
BY  SELECTING  SENSOR  LOCATIONS  TO  MAXIMIZE  U. 


- NEXT  SLIDE  PLEASE  - 
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For  low  output  i 
V((o) 

For  high  output 
V(w) 


>NR 


= M^  S^(gj) 
n2(cj) 


SNR 


= M S(cu) 
N(o;) 
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The  denominator  term^  V of  omega^  is  a complicated  function 

OF  SIGNAL-TO-NOISE  RATIO  (SNR)  AND  THE  NUMBER  OF  SENSORS^ 

The  output  signal-to-noise  ratio  is  defined  by  the  product  of 

THE  INPUT  SIGNAL-TO-NOISE  RATIOy  S OF  OMEGA  OVER  N OF  OMEGA^ 
multiplied  by  the  NUMBER  OF  SENSORS^  ThE  FUNCTION^  V OF 
OMEGA^  BEHAVES  DIFFERENTLY  AT  LOW  AND  HIGH  OUTPUT  SNR,  FOR  LOW 
OUTPUT  SNR^  THIS  FUNCTION  BEHAVES  AS  THE  SQUARE  OF  THE  OUTPUT 

SNR.  Assuming  we  are  in  a passive  situation  and  cannot  con  ^ol 

THE  INPUT  SIGNAL-TO-NOISE  RATIO^  WE  SEE  OUR  ONLY  METHOD  FOR 
IMPROVEMENT  IN  THE  V FUNCTION  IS  TO  INCREASE  M.  WhEN  WE  ARE  AT 
LOW  OUTPUT  SNR,  INCREASING  M HELPS  DRAMATICALLY.  HoWEVER,  AS 
M IS  INCREASED,  THE  OUTPUT  SIGNAL-TO-NOISE  RATIO  GETS  LARGER  AND 
THE  V FUNCTION  IMPROVES  LESS  DRAMATICALLY;  INDEED,  INCREASING 
M IS  LIKE  INCREASING  THE  INTEGRATION  TIME. 


- NEXT  SLIDE  PLEASE  - 
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The  term  U that  we  desire  to  maximize  is  different  for 

RANGE  AND  BEARING  ESTIMATION.  FOR  BEARING  ESTIMATION  WITH  FOUR 
PROPOSED  LINE  ARRAYS^  U SUB  B IS  GIVEN  BY  THE  SQUARE  OF  THE 
EFFECTIVE  ARRAY  LENGTH  DIVIDED  BY  A CONSTANT.  ThE  EFFECTIVE 
ARRAY  LENGTH  IS  EQUAL  TO  THE  PHYSICAL  ARRAY  LENGTH  TIMES  THE  SINE 
OF  THE  BEARING  TO  THE  SOURCE.  ThUS^  IT  IS  THE  APPARENT  LENGTH  OF 
THE  ARRAY  AS  SEEN  BY  THE  SOURCE.  In  BEARING  ESTIMATION^  WE  DE- 
SIRE TO  MAKE  L SUB  E LARGE  AND  THE  CONSTANT  K SUB  B SMALL  IN 
ORDER  TO  REDUCE  VARIANCE.  NoTE  THAT  DOUBLING  THE  ARRAY  LENGTH 
REDUCES  THE  VARIANCE  BY  FOUR.  ThUS^  ARRAY  LENGTH  IS  A MORE  IM- 
PORTANT FACTOR  IN  BEARING  ESTIMATION  THAN  EITHER  INTEGRATION  TIME 
OR  THE  NUMBER  OF  HYDROPHONES  WHEN  ONE  IS  ALREADY  AT  HIGH  OUTPUT 
SIGNAL-TO-NOISE  RATIO. 

The  four  different  array  types  studied  are  an  equispaced 

LINE  ARRAY  AND  THREE  ARRAYS  WITH  M ELEMENTS  GROUPED  AT  THE  TWO 
ENDS  AND  the  MIDDLE  OF  THE  ARRAY.  McDoNALD  AND  SCHULTHEISS  HAVE 
SHOWN  THAT  BY  PLACING  HALF  OF  THE  M ELEMENTS  AT  EACH  END  OF  A 
LINE  ARRAY  IN  AN  FI  OVER  TWO^  ZERO^  FI  OVER  TWO  GROUPING^  A BOUND 
ON  BEARING  VARIANCE  IS  OBTAINED.  ThREE  GROUPS  OF  A THIRD  EACH 
RESULT  IN  A K SUB  B OF  THREE  AND  GROUPINGS  OF  QUARTER^  HALF, 
QUARTER  RESULT  IN  A K SUB  B OF  FOUR. 

The  practical  implications  of  FIcDonald  and  Schultheiss's 

RESULT  ARE  BOTH  TO  PROVIDE  A BOUND  ON  HOW  WELL  BEARING  COULD  BE 
ESTIMATED  UNDER  IDEAL  CONDITIONS  AND  TO  SUGGEST  HOW  TO  PLACE  A 
LIMITED  NUMBER  OF  HYDROPHONES  OVER  A LARGE  APERTURE.  NOTABLY 
HALF  OF  THE  HYDROPHONES  SHOULD  BE  POSITIONED  AT  EACH  END  OF  THE 
ARRAY,  PLACED  AT  HALF-WAVL  LENGTH  SPACING  FOR  THE  DESIGN 
FREQUENCY. 


-NEXT  SLIDE  PLEASE- 
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The  term  U for  ranging  depends  on  the  fourth  power  of  the 

RANGE  RELATIVE  TO  THE  EFFECTIVE  BASELINE.  ThE  VARIANCE  OF  THE 
RANGE  ESTIMATE  IS  REDUCED  BY  MAKING  THE  EFFECTIVE  ARRAY  LENGTH  L 
SUB  E LARGE.  ThIS  CAN  BE  DONE  BY  MAKING  THE  ARRAY  LENGTH  L 
LARGE  OR  BY  PHYSICALLY  STEERING  THE  ARRAY  BROADSIDE  TO  THE  SOURCE. 

The  variance  can  also  be  reduced  by  decreasing  the  range  to  the 
SOURCE.  Most  notable  in  the  U sub  R term  is  that  the  variance  of 

THE  RANGE  ESTIMATE  INCREASES  AS  THE  FOURTH  POWER  OF  THE  RATIO  OF 
RANGE  TO  THE  EFFECTIVE  ARRAY  LENGTH.  ThIS  FOURTH  POWER  DEPENDENCE 
OF  THE  VARIANCE  ON  ARRAY  LENGTH  SUGGESTS  THE  NEED  TO  MAKE  L SUB  E 
LARGE  IN  ORDER  TO  ESTIMATE  RANGE. 

The  constant  K sub  R depends  on  the  array  type.  For  an 
EQUISPACED  LINE  ARRAY^  K SUB  R IS  360.  A BOUND  IS  PROVIDED  BY 
AN  ARRAY  CONFIGURED  WITH  A QUARTER  OF  THE  HYDROPHONES  AT  EACH 
END  AND  HALF  IN  THE  MIDDLE.  ThUS,  WE  SEE  THAT  THE  ARRAY  CON- 
FIGURATION DESIRED  FOR  BEARING  ESTIMATION  AND  THE  ONE  FOR  RANGE 
ESTIMATION  DIFFER.  ThE  BEARING  ARRAY  SHOULD  HAVE  ITS  ELEMENTS 
TOWARD  THE  ARRAY  ENDS^  WHILE  THE  RANGING  ARRAY  SHOULD  HAVE  HALF 
OF  ITS  ELEMENTS  IN  THE  CENTRAL  PORTION. 

As  POINTED  OUT  IN  THE  CONFERENCE  PROCEEDINGS,  AN  ARRAY  WITH 
A THIRD  0^ ITS  ELEMENTS  AT  EACH  END  AND  THE  MIDDLE  WILL  MINIMIZE 

Thus,  in  some  sense  an  array  comprised  of 

THREE  EQUAL  GROUPS  OF  ELEMENTS  WILL  OUTPERFORM  ALL  OTHER  ARRAYS 
FOR  PASSIVELY  LOCATING  AN  ACOUSTIC  SOURCE. 

- NEXT  SLIDE  PLEASE  - 
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The  relative  range  error  given  by  the  standard  deviation  of 

THE  RANGE  ESTIMATE  DIVIDED  BY  THE  TRUE  RANGE  IS  GIVEN  BY  A CONSTANT 
TIMES  THE  STANDARD  DEVIATION  OF  THE  BEARING  ESTIMATE  (MEASURED 
IN  radians)  times  a term  that  DEPENDS  LINEARLY  ON  THE  RANGE  TO 
THE  SOURCE  RELATIVE  TO  THE  EFFECTIVE  ARRAY  LENGTH.  It  IS  IMPORTANT 
TO  NOTE  THAT  THE  RELATIVE  ERROR  INCREASES  LINEARLY  WITH  ACTUAL 
RANGE  TO  THE  SOURCE.  STANDARD  DEVIATIONS  OF  RANGE  ESTIMATES 
INCREASE  AS  THE  SQUARE  OF  THE  RANGE  TO  THE  SOURCE.  HeNCE  WE  SEE 
THAT  IT  IS  EXTREMELY  DIFFICULT  TO  ESTIMATE  RANGE  OF  A DISTANT 
SOURCE  EVEN  UNDER  IDEAL  CONDITIONS. 

One  of  the  advantages  of  expressing  relative  range  errors 

IN  THIS  FORM  IS  THAT  THE  STANDARD  DEVIATION  OF  BEARING  ESTIMATES 
IS  A TERM  FAMILIAR  TO  SONAR  ENGINEERS  AND  SIGNAL  PROCESSORS. 

Moreover^  there  is  concern  that  the  ocean  medium  may  inherently 

LIMIT  THE  PRACTICAL  ABILITY  TO  ESTIMATE  BEARING  EVEN  THOUGH  THEORY 
PREDICTS  WITH  ENOUGH  SIGNAL-TO-NOISE  RATIO  OR  INTEGRATION  TIME^ 

THE  BEARING  CAN  BE  MEASURED  ARBITRARILY  WELL.  ThE  EXPRESSION 
GIVEN  HERE  CLEARLY  POINTS  OUT  THE  NEED  TO  MAKE  THE  ARRAY  LENGTH 
LARGE  WHEN  THE  SOURCE  RANGE  CANNOT  BE  REDUCED.  Of  INTEREST  IS 
THAT  THIS  CONCLUSION  IS  EXTREMELY  INSENSITIVE  TO  THE  TYPE  OF 
ARRAYy  PROVIDED  THE  ARRAY  HAS  SOME  RANGING  CAPABILITY.  ThIS 
CAN  BE  SEEN  FROM  THE  SIMILARITY  OF  THE  CONSTANTS  GIVEN. 


- NEXT  SLIDE  PLEASE  - 
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To  SUMMARIZE,  WE  DESIRE  TO  KNOW  HOW  TO  PLACE  A LIMITED 
NUMBER  OF  HYDROPHONES  OVER  A SPECIFIED  BASELINE.  ThE  HYDRO- 
PHONES SHOULD  BE  PLACED  IN  GROUPS,  WITH  THE  HYDROPHONES  IN  EACH 
GROUP  PLACED  AT  HALF  WAVE  LENGTH  SPACING  FOR  THE  DESIGN  FREQUENCY. 

For  bearing  estimation,  half  of  the  FI  hydrophones  should  be 

PLACED  AT  EACH  END  OF  THE  ARRAY.  FOR  RANGE  ESTIMATION,  A QUARTER 
OF  THE  HYDROPHONES  SHOULD  BE  PLACED  AT  EACH  END  OF  THE  ARRAY  AND 
HALF  PLACED  IN  THE  MIDDLE.  FOR  SIMULTANEOUSLY  ESTIMATING  RANGE 
AND  BEARING,  THE  HYDROPHONES  ARE  PLACED  IN  THREE  GROUPS  EACH 
WITH  M OVER  THREE  HYDROPHONES. 

SLIDE  OFF  PLEASE.  ARE  THERE  ANY  QUESTIONS? 
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MfTHODS  FOR  PASSIVELY  LOCATING  AN  ACOUSTIC  SOURCE 


Or.  G.  Clifford  Carter 


Naval  Underwater  Systems  Center.  New  London,  CT  06320 


ABSTRACT 

The  maximum  likelihood  (ML)  processor  is  pre- 
sented for  passively  estimating  range  and  bearing 
to  an  acoustic  source.  The  source  signal  is  ob- 
served for  a finite  time  duration  at  several  sen- 
sors in  the  presence  of  uncorrelated  noise  When 
the  speed  of  sound  in  an  isovelocily  medium  and 
the  sensor  positions  are  known  the  ML  estimator 
for  position  constrains  the  source  to  sensor  de- 
lays to  be  focused  into  a point  corresponding  to 
a hypothesised  source  location.  The  variances  of 
the  range  error  and  bearing  error  are  presented 
for  the  optimum  processor.  It  is  shown  that  for 
bearing  and  range  estimation  different  sensor 
configurations  are  desirable.  However,  if  the 
area  of  uncertainty  is  to  be  minimized,  then  the 
sensors  should  be  divided  into  equal  groups  with 
one  third  of  the  sensors  in  each  group. 


INTRODUCTION 

An  (underwater)  acoustic  point  source  radiat- 
ing energy  to  four  (or  more)  receiving  sensors 
is  shown  in  Figure  (1).  The  position  of  the 
source  (in  two  space)  can  be  characterized  by  po- 
lar coordinates  (range  and  bearing)  for  a given 
frame  of  reference. 


kr.;. 


length  I ’ 

FIGURE  (1)  R,  8.  L.  M Defined 

For  example,  the  range,  R,  and  bearing,  B.  rela- 
tive to  a specified  sensor  and  line  from  the  sen- 
sor can  be  used  as  shown  in  figure  (1).  The 


problem  we  address  here  is  how  to  estimate  range, 
R,  and  bearing,  B,to  a source  when  M sensors 
separated  by  a maximum  of  f (meters)  have  observed 
T seconds  of  received  data.  Ue  will  examine  the 
maximum  likelihood  (ML)  technique  for  position 
estimation. 

MATHEMATICAL  MODEL 

For  our  purposes  we  assume  that  each  receiving 
sensor  at  the  t-th  instant  in  time  corresponds  to 
a signal  plus  noise.  Namely,  the  i-th  (of  M)  sen- 
sor outputs  is  characterized  by 

(1) 


rj(t)  = 
0<t<T. 


Si  (t)  + ni  (t ) . i = 1 , M 


The  signal  and  noise  are  uncorrelated  and  the 
noises  are  mutually  uncorrelated. 

For  a spatially  stationary  source  the  signal 
can  be  viewed  as  a (range  and  frequency  dependent 
attenuated  and)  delayed  source  signal.  More  gen- 
erally, source  motion  will  impart  a time  compres- 
sion (or  expansion)  to  the  received  signal  wave 
foiTii.  Source  motion  is  discussed  by  Knapp,  Carter 
and  others,  in  this  conference,  but  is  beyond  the 
scope  of  this  paper.  However,  we  know  ttie  prob- 
lem of  estimating  the  position  of  a stationary 
source  is  considerably  easier  than  that  of  a mov- 
ing source.  Thus,  the  lesults  in  this  paper  serve 
as  a bound  on  performance  (variance),  still  we 
will  see  that  it  is  extremely  difficult  in  the 
best  case  to  estimate  source  range.  We  inow  fur- 
ther (for  a moving  source)  if  the  time  compression 
varies  (significantly)  across  the  sensor  baseline 
that  failure  to  compensate  for  source  motion  will 
(significantly)  degrade  performance.  Extracting 
estiniates  of  time  compression  may  aid  in  motion 
analysis  and  indeed  may  he  necessary  but  certain- 
ly performance  would  be  l>etter  against  a station- 
ary source  (or  moving  source  with  known  motion). 

While  in  general,  attenuation  is  a function  of 
frequency  and  range,  we  can  initially  treat  atten- 
uation as  a constant  and  later  impose  the  frequen- 
cy (and  range)  dependence.  In  particular,  we  can 
model  the  i-th  signal  as 

sj(t)  = ajS(t  ♦ D,),  i = 1,  M (?) 
where  s(t)  is  the  signal  at  the  source. 

(Without  loss  of  generality,  we  can  assume  ’ 1 
and  D,  - 0. ) 

There  has  been  considerable  attention  devoted 
to  estimating  the  (relative)  time  delay  parameter 
vector  (D«  - Di , Dj  - Di , . , D|,^  - Ui ) ^ee  for 

example,  MacDonald  and  Schultheiss  (lybd).  Bangs 
and  Schultheiss  (1973),  Hannan  and  Thompson  0 573), 
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ridhn  (1975). anc  »’r,  ..p  and  Carter  {1976). 
c,AF:‘(G  ISTIMATION 

The  relative  delay  between  two  receiver  (or 
bpdinfo mpr ) outputs  is  related  to  the  bedring  of 
the  source.  In  particular,  the  bearing  of  the 
source  relative  to  ‘he  i - th  and  (i  - 1)-St 
sensor  is 


where  C is  the  speed  of  sound  in  a nondi spersive 
isovelocity  medium  and  d is  the  separation  be- 
tween the  two  sensors  (or  beainformer  centers). 

The  bearing  relative  to  the  defined  reference  line 
is  known  when  the  bearing  to  the  source  relative 
to  the  i - th  and  (i  - l)st  sensor  is  known  and 
the  position  of  those  sensors  is  known  (relative 
to  the  reference  line).  In  particular,  the  effec- 
tive bearing  angle  is  v/hat  you  would  see  at  two 

sensors  corrected  by  how  the  two  sensors  are  out 
of  alignment  (relative  to  the  reference  line).  In 
practice,  errors  in  bearing  are  caused  by  estima- 
tion errors  and  by  uncerta i ni t i es  in  knowing  the 
sensor  location.  If  the  sensor  positions  vary 
with  time  about  some  known  position,  time  averag- 
ing will  correct  for  position  variation.  If  the 
sensor  positions  vary  about  some  unknown  position, 
significant  biases  in  bearing  (range)  estimation 
can  exist  that  cannot  be  averaged  out.  The  esti- 
mation errors  are  a function  of  integration  time, 
processing  bandwidth  and  coherence  (or  signal -to- 
ne ise  ratio  (SNR))  between  the  two  sensors,  as 
well  as  inhomogen i t i es  in  the  medium  causing  wave 
front  warping  (difference  in  C between  the  source 
and  various  sensors).  It  can  be  shown  that  (for 
many  practical  problems)  we  can  transform  our 
pioblrm  into  an  cguivalent  one  in  wnich  the  range 
and  bearing  to  a source  are  to  be  measured  by 
measuring  the  effect ’ve  bearings  to  a linear  array. 
Thus,  we  concentrate  our  attention  on  linear  ar- 
rays. When  the  cfler’.ive  be. w ings  ai^e_k;^nowr^  c>^- 
a_c  t ly . four  or  more  sensor  bearing  Tpairs)  provide 
a redundant  characlerizat ion  of  the  source  posi- 
tion. If,  however,  we  estim.ite  the  bearings  and 
plot  the  bearing  lines  (or  asyiiptoles  to  hyper- 
bolic curves)  then,  in  general,  the  lines  will  no_t 
intersect  in  a single  point. 

From  the  discussion  of  thp  mathematical  model, 
we  can  see  a method  evolving  to  estimate  the  posi- 
tion of  the  source.  This  method  focusi.s  all  bear- 
ing lines  into  one  hypothesited  source  position. 

KAXIMUM  I IKI 1 I HOOD 

The  maximum  likelihood  ('ll)  estimate  for  the 
lime  delay  vector  has  been  derived  fur  stationary 
Gaussian  processes  by  Hahn  ( 1975)  and  Carter  ( 1976) . 
As  stated  earlier,  these  time  delay  estimates  give 
rise  to  several  (M-1)  hearing  lines  (ft  cm  (3)) 
which  do  not  necessarily  int'rsect  in  one  point. 

It  is  conceivable  that  a point  that  is  the  (weight- 
ed) least  squares  distance  to  all  bearing  lines  is 
mt  the  best  {maximum  likelihnod)  estimate  of 
source  position.  It  is  not  difficult  to  show  that 
ML  estimate  for  range  and  bearing  when  the  sensor 
(element)  positions  arc  known  is  achieved  by  a 


simple  variation  on  Hahn  (1975)  and  Carter  (1976). 
In  particular,  by  focusing  all  the  time  delay  ele- 
ments at  many  (hypothesized)  range  and  bearing 
pairs  and  watching  for  the  peak  output  of  the  ML 
time  delay  vector  system,  the  ML  position  estimate 
is  observed.  Another  way  to  look  at  this  problem 
is  that  we  want  to  maximize  a quantity  by  adjust- 
ing a number  of  delay  parameters  subject  to  the 
constraint  that  all  the  delays  must  focus  (inter- 
sect) in  a (nominally)  single  (hypothesized)  posi- 
tion. Equivalently,  rather  than  adjust  M-1  delay 
dials  and  then  have  to  figure  out  a least  squares 
position  from  them,  we  can  adjust  one  range  and 
one  bearing  dial,  and  by  a simple  (perhaps  pre- 
computed) trigonometric  mapping,  select  all  (M-1) 
hypothesized  delays  necessary  to  focus  the  sensors 
at  one  point  in  two'  space.  The  performance  of 
Such  a system  has  been  examined  by  Bangs  and 
SchuUheiss  (1973). 

For  constant  attenuation,  when  (and  only  wtien) 
we  utilized  the  (proper  frequency  dependent)  pre- 
filters given  by  (see  Bangs  and  SchuUheiss  ( 1973)), 
Hahn  (1975),  Carter  (1976),  and  Knapp  and  Carter 
(1976). 

S(m,^)/N2(m|,) 

(where  S and  N denote  signal  and  noise  power  spec- 
tral densitie;),  then  the  variance  of  the  parameter 
6 (where  G=R  orO=B)ls  ^ 

(0)»nl|ETr(P0P0)T/^u^:S(w^.)h|(U,^]  (5) 

where  * denotes  complex  conjugate  transpose,  Tr 
denotes  the  trace,  P.  is  a particular  position 
matrix,  and  E is  an  effective  degradation  factor 
both  to  be  defined,  for  bearing  estimation,  the 
i - th,  j - th  element  of  the  Pg  matrix  is  given 
by 

= 5Jn  B (2-  . 2.)  (6) 


where  Z 

sensor,  i or  range  estimation,  the  i 
element  of  the  Pj,  matrix  is  given  bv 

ii?  - 


- th 
- th,  j - th 


Pij 


is  defined 


(7) 


The  degradation  factor 

Tr(PBpB)1r(PRPR)  -[ir(pBPR)] 


(8) 


E = 


lr(PBPg)lr(PgPf() 


Thus,  when  (but  only  when)  7r(PgPg)=  0 we  see  that 
E = 1 and  there  is  no  degradation . However,  in 
general,  this  will  not  be  the  case  and  E must  be 
computed.  Pi oper  (or  forlui tous )choice  of  the 
coordinate  system  will  insure  E ^ I and  then  only 
Tr(P„P!)  plays  a role  in  the  expression  for  vari 


'00 

ance  of  the  estimates.' 


An  important  te 

S(ujhk 


n the  variance  expression  is 

_aJL  (5) 

1 ^ tta 
where  a - S{ui|j )/N(u|^ ) 


At  first  glance,  it  may  appear  that  the  variance 


♦Errata : We  can  show  for  a synmetric  line  array  (SLA),  with  the  origin  in  the  center, 

that  E=l.  All  the  results  that  follow  apply  to  bearings  and  ranges  measured  to 
the  center  of  an  SLA.  A more  complete  discussion  will  be  given  in  a subsequent 
paper. 
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could  be  reduced  by  increasing  h;  however,  if  the 
sensor  positions  are  fixed  then  the  variance,  giv- 
en by  (5),  is  the  minimum  variance  obtainable  and 
is  attained  when  h)(  is  of  tiie  form  specified  by 
(4).  For  low  output  SNR  (Ma  < < 1)  (9)  goes  as 
the  square  of  the  SNR.  For  high  output  SNR  (9) 
goes  as  the  ratio  of  SNR  to  M. 

When  the  number  of  sensors  and  the  SNR  fix  the 
output  SNR,  the  performance  (variance)  is  governed 
only  by  CTTr(P.  P*).  Increasing  the  observation 
time,  T,  decreases  the  variance  according  to  (5). 
We  desire  to  know  how  to  decrease  the  variance 
through  proper  element  placement.  MacDonald  and 
Schultheiss  (1969)  have  stated  that  for  bearing 
estimation  half  of  the  elements  should  be  placed 
at  each  end  of  an  array  of  fixed  length  (provided 
the  noise  remains  uncorrelated).  In  practice,  if 
we  are  constrained  to  have  M elements  in  an  array 
of  length  L,  we  should  place  M/2  by  beginning  at 
each  end  and  spacing  every  half  wave  length  (for 
the  design  frequency)  towards  the  center. 

For  an  eqiiispaced  line  array  (d  spacing)  we 
find  for  large  M and  L ^ (M-l)d 

[ETr(PgP;)]  = 6cL  (,o) 

m2 


range  or  bearing  is  to  be  estimated).  If  we  de- 
fine that  the  average  error  vector  be  minimized 
between  our  position  estimate  and  the  true  source 
position  then  it  can  be  shown  that  the  sensor 
positions  will  be  a function  of  the  true  range  to 
the  source.  (In  particular,  for  close-in  sources 
we  will  want  a (M/2,0,  M/2)  division  of  sensors 
and  for  very  distant  sources  we  will  want  element 
groups  of  M/4,  H/2,  M/4.)  One  way  of  defining  a 
best  array  geometry  is  to  minimize  the  area  of  un- 
certainty. That  is,  select  the  sensor  positions 
to  minimize  the  quantity  Ro(B)a(R).  Such  an  area 
of  uncertainty  will  clearly  not  be  square;  for 
example,  at  large  ranges  we  will  bo  more  apt  to 
make  range  errors  than  bearing  errors.  The  opti- 
mum placement  was  found  numerically  to  be  M/3, 
M/3,  M/3.  For  the  one  third,  one  third,  one  third 
division 


ni) 

m2  (14)* 


where  the  effective  array  length  projec- 
tion seen  by  the  source  is  given  by  Le  = LsinB, 
which  agrees  with  eq(34)  of  Gangs  and  Schultheiss 
(1973)  and  agrees  with  eq(19)  of  MacDonald  and 
Schultheiss  (1969).  As  expected,  good  perfor- 
mance Oow  variance)is  obtained  by  steering  broad- 
side to  the  source  and  making  the  array  length.  L, 
and  number  of  sensors,  M,  as  large  as  possible. 

In  theory,  the  variance  can  be  reduced  by  a factor 
of  three  by  proper  element  placement  (That  is.  the 
6 in  (10)  is  replaced  by  a 2). 

The  ranging  perfonnance,  for  an  equispaced 
line  array  can  he  shown  to  be  governed  by 


3fc0^  C3p4 

t V 


(11)  * 


This  result  agrees  with  Gangs  and  Schultheiss 
(1973).  In  expression  (11)  we  see  the  tremendous 
importance  of  absolute  range  to  the  source  (vari- 
ance goes  as  a fourth  power).,,  In  order  to  corrob- 
orate the  derivation,  Flr(P„Pp)  was  numerically 
computed  and  the  equispaced  results  ac|reed  with 
the  derivation  as  expected.  Nuttall  (1976)  has 
shown  that  if  the  array  elements  were  symmetrical- 
ly placed  about  scx-ie  origin  that  the  best  ranging 
perfornmnee  (lowest  variance)  could  be  achieved  by 
putting  one  half  the  elements  in  the  middle  and 
one  quarter  at  each  end  of  tie  array.  (However, 
there  can  bo  non  synr.ietric  ranging  arrays  that  out- 
perform the  sycinetric  array.)  For  an  array  of 
quarter,  half,  quarter  grouping  it  can  be  shown 
that 


E'--(PrPR>] 


-1 


IZSc-rI 

iy~ 


(12) 


which  outperforms  the  equispaced  array  by  a factor 
of  2.8  in  variance,  and  is  conjectured  to  be  the 
best  ranging  array  geometry  for  large  M. 

There  is  a difference  between  how  to  place  M 
elements  of  a linear  array  depending  on  v/hether 


Thus,  by  placing  M/3  sensor  elements  at  -L/2,  M/3 
elements  at  0 and  M/3  elements  at  L/2  we  achieve 
a comprise  bearing  and  range  estimation  system, 
such  a comprise  minimizes  the  area  of  uncertainty 
defined  by  Ro(R)o(B).  Attempts  to  minimize  other 
error  measures  can  result  in  array  configurations 
that  depend  on  source  range. 
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♦Errata;  The  coefficients,  360,  128  and  144,  in  equations 
here  are  believed  to  be  correct;  they  are  16  times  larger 
original  manuscript. 
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than  reported  in  the 
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ERRATA  ANO  APPENDIX 

METHODS  FOR  PASSIVELY  LOCATING  AN  ACOUSTIC  SOURCE 
Dr.  G.  Cl ifford  Carter 


In  practice,  if  we  had  a limited  number  of  elements, 
they  snould  be  placed  at  half  wavelength  spacing 
(for  a given  design  frequency)  in  three  groups, 
with  M/3  at  each  end  of  the  baseline  and  M/3  in 
the  center.  It  is  useful  to  note  for  three  groups 
of  M/3  that  theoretically 


1.  Comparing  eq.  (4)  and  (9)  we  see  an  inconsist-  g(R)  . , * R 

ency.  Equation  (4)  should  be  R ' ° 


h 


k 


S(ui|^)/N^(.a,^) 


where  h|(  is  the  magnitude  squared  transfer  function. 


where  c(B)  is  measured  in  radians.  Increases  in 
M,  T or  SNR  decrease  o(B)  and  thus  only  indirectly 
decrease  the  relative  range  error.  Increases  in 
Lg  effect  both  o(B)  and  the  relative  range  error. 
U is  also  instructive  to  look  at  the  ratio  of  the 
major  and  minor  axes  of  the  uncertainty  ellipse. 

We  note 


2.  If  in  eq.  .5)  we  define 


2~t 

T 


= 4 /J  ?! 
o(B)  e 


as  is  usually  done,  then  equation  (5)  must  be  Thus,  the  shape  of  the  uncertainty  ellipse 

multiplied  by  2-  in  order  to  Oe  correct.  The  changes  depending  on  range  and  effecti'.e  baseline. 

rationale  fur  tne  form  of  (5;,  which  is  mcd’fied 

from  Bangs  and  Schulrheiss  (1973),  is,  of  course, 

to:  (1)  show  the  dependence  on  integration  time 

and  (?)  allow  tre  sutiwanon  to  be  replaced  by 

integral  cn.  It  should  also  be  oointed  out  that 

eq.  (5);  holds  only  for  high  output  SNR,  involves 

the  local  variat'on  of  the  parameter  estimate  and, 

ignores  a .biguous  paramete'-  estimates  due,  for 

examp’e.  to  ’arge  sidelobes. 

3.  Following  the  discussion  (about  E-1  near 
eq.  (9))1  a comment  should  have  tee.n  inserted 

that  "we  can  show  for  a s.nnetric  I'r.e  array  (SI A), 
with  the  origin  in  the  ce'^te-',  that  E-1.  All  the 
resulrs  that  follow  apply  to  bearing  ano  ranges 
measured  to  the  center  o'  an  SLA  and  not  to 
Figure  1 of  the  text.  A nore  como'ete  discussian 
will  be  given  in  a subseouert  pater." 

4.  the  coefficients  in  equations  ;il),  (12)  and 
(14)  snould  be  360,  123  and  144  respe.rtively; 
tnat  is,  the  ranging  vari.<rre  is  16  times  worse 
than  queted  i.n  ‘he  uncorrected  manussrict. 

5.  The  degradation  ‘acor,  eg.  (S),  can  be 

viewed  as  follows.  For  any  array  ard  any  source 
a fixed  el ’ipse  of  uncertainty  ("ists  fer  thg 
particular  fixed  SNR,  L.  S,  and  R.  The 

absolute  position  and  area  of  this  ellipse  da  not 
change  as  a function  of  :ne  coordinate  system. 

However,  if  when  we  select  tne  coordinate  system 
we  define  me  range  and  bearing  as  being  measured 
from  the  center  o'"  'ne  coordinate  system  then  the 
variance  of  t^e  range  and  bearing  estimates  do 
cnange  .js  a function  o'  the  coordinate  system. 

6.  Tne  minimum  area  if  uncertainty  is  not 
necessarily  obMinable  thus  equations  (13^  and 
(14)  serve  as  a bouno  on  perfornance, 

and  further,  do  not  consider  issues  o'"  beam 
pattern  sidelobes  causing  ambiguous  estimates. 
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